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In a recent paper Hao et al. [Phys. Rev. B 80, 035329 (2009)] reported variational calculations 
of energy spectrum for shallow hydrogenic donor in the structure of semiconductor/insulator/metal 
with a new type of trial wave function. They also performed calculations for semiconductor/insulator 
system and found that their method gives energy values lower than those obtained by MacMillen and 
Landman [Phys. Rev. B 29, 4524 (1984)]. As follows from these results MacMillen and Landman 
have got much larger errors in energy values than they expected. However we confirm that the 
theoretical approach suggested by MacMillen and Landman gives rather accurate energy values for 
the system of hydrogenic donor near the interface between semiconductor and insulator. 



In a recent paper—, Hao et al. (HDAP) reported 
variational calculations of energy spectrum for shal- 
low hydrogenic donor in the structure of semiconduc- 
tor/insulator/metal with a new type of trial wave func- 
tion. To evaluate the quality of the suggested trial func- 
tion HDAP compared the calculated ground state en- 
ergies at different donor positions with the calculations 
obtained by MacMillen and Landman^. The authors of 
Ref. @ ensured that their calculated energy values are 
accurate to four significant figures. However, in Ref. [l[ 
it was claimed that the errors in the result of Ref. @ may 
be even grater than 9%. The aim of this comment is to 
resolve this contradiction. We argue that the comparison 
in Ref. [l| is not adequate. 

To calculate the electronic spectrum under consider- 
ation an eigenvalue problem for the following equation 
(written in dimensionless form) should be solved 2 -: 
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F(r) = EF(f). (1) 



All distances are scaled to units of effective Bohr radius 
(a* B ) and the energy E is given in effective Rydbergs 
(Ry*). The parameter Q* is given by: 
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where e s and Ci ns are the dielectric constants of the semi- 
conductor and the insulator respectively. The distance of 
the donor from the interface is R and and tq are the 
distances of an impurity electron from the donor (D) and 
its image (D), respectively. 

HDAP compared their results calculated for Q* = 1 in 
Eq. ([1]) (see Table I in Ref. [1[) to the values calculated 
in Ref. for Q* w 0.8387 (e s = 11.4, e ms = 1). The 
latter one corresponds to the interface between silicon 
and vacuum (see Table II in Ref. 2]). It seems that this 
incorrectness has been arisen due to the same mistake 
made in Ref. Q. 



In order to test this assertion we have performed cal- 
culations for Q* « 0.8387 using variational and finite- 
element approaches. For variational calculations a trial 
wave function was chosen as a sum Atrial = ^Cik&k 
with the basis wave functions: 



= 2a 3 / 2 i?^ 1 / 2 exp [-aR(t - r,)} Li{t,)P k (rj), (2) 



where £ = 



rn+ri 



and rj = 



are prolate spheroidal 
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coordinates, a is a variational parameter, and and 
Pk(v) are polynomials of i-th and fc-th degrees, respec- 
tively. As Pk{f]) we choose the Legendre polynomials of 
odd degree which allow to satisfy the boundary condition 
F(r) = at the interface (z = 0). Choosing this basis 
we avoid the necessity of numerical integration. Another 
advantage of this basis set is a lower number of basis 
functions which is necessary to obtain the same accuracy 
for the ground state energy Eq as in Ref. @. 

To carry out numerical computations using finite ele- 
ment method (FEM), the system is assumed to be rota- 
tionally symmetric around z-axis and the problem is re- 
formulated in cylindrical coordinates (r, z). The equation 
is posed in a bounded domain Q, = (0, <5) x (—8 — R, 0) 
with S = 10 chosen for an approximation of the semi- 
infinite region (0, +oo) x (— oo,0) with a donor position 
at (0, — R). Using MATLAB we discretize this problem 
by linear finite elements on a triangular mesh with a 
number of unknowns w 150000. The generalized sparse 
eigenvalue problem is solved by the implicitly restarted 
Arnoldi method in MATLAB. 

The comparison of our results for Q* « 0.8387 with 
the theoretical results of MacMillen and Landman is pre- 
sented in Table U which has the form similar to that one 
in Ref. [l[ . We found that our values are practically co- 
inciding to those obtained by MacMillen and Landman 
at the same distances from the interface. So our calcula- 
tions confirm the evaluation of the significant figures for 
energy values presented by MacMillen and Landman. 

To evaluate the quality of the trial function with two 
parameters suggested by HDAP we have also performed 
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TABLE I. The values of the ground state energy Eo of the 
impurity electron near the semiconductor-insulator interface 
for different values of R/a* B calculated on the basis of the 
variational wave functions of Eq. (J2j) with Q* — 0.8387, and 
Eo found in Ref. @] for the same value of Q* . 



R/a* B 


a 


E /Ry* 


E /Ry* 


Eo/Ry* 






variational, 


FEM, 


variational, 






present 


present 


Ref. [2J 


0.2 


0.854 


-0.6064 


-0.6062 


-0.6064 


0.4 


0.908 


-0.6508 


-0.6507 


-0.6507 


0.6 


0.961 


-0.7223 


-0.7222 


-0.7221 


0.8 


0.988 


-0.8099 


-0.8098 


-0.8098 


f.O 


0.908 


-0.8946 


-0.8946 


-0.8945 


f .2 


0.827 


-0.9640 


-0.9643 


-0.9640 


1.4 


0.773 


-1.0158 


-1.0164 


-1.0158 


1.6 


0.639 


-1.0522 


-1.0530 


-1.0521 


1.8 


0.666 


-1.0767 


-1.0771 


-1.0767 


2.0 


0.666 


-1.0925 


-1.0932 


-1.0925 


3.0 


0.666 


-1.1085 


-1.1089 


-1.1086 


4.0 


0.666 


-1.0943 


-1.0952 


-1.0944 


5.0 


0.720 


-1.0794 


-1.0804 


-1.0794 


6.0 


0.720 


-1.0676 


-1.0679 


-1.0676 



TABLE II. The values of the ground state energy Eo of the 
impurity electron near the semiconductor-insulator interface 
for different values of R/a* B calculated on the basis of the 
variational wave functions of Eq. ([2]) with Q* = 1, and Eo 
found in Ref. [l( for the same value of Q* . 



R/a* B 


Eo/Ry* 


Eo/Ry* 


Eo/Ry* 


Realtive 




variational, 


FEM, 


variational, 


error 




present 


present 


Ref. [1] 


in% 


0.4 


-0.7208 


-0.7208 


-0.716 


0.6 


1.0 


-0.9491 


-0.9491 


-0.927 


2.3 


1.6 


-1.0917 


-1.0925 


-1.077 


1.4 


2.0 


-1.1256 


-1.1264 


-1.116 


0.9 


3.0 


-1.1323 


-1.1327 


-1.128 


0.4 


4.0 


-1.1131 


-1.1139 


-1.111 


0.3 


6.0 


-1.0806 


-1.0810 


-1.080 


< 0.1 



variational and FEM calculations for Q* = 1. The ob- 
tained results are presented in Table HU As seen from this 
table the trial function suggested by HDAP gives rather 
good bound energy values. However, they are less than 
the values obtained using the trial wave function consist- 
ing of a sum over a basis set of wave functions (as in Ref. 
2] or ours) or by FEM. 

HDAP also examined the quality of their trial func- 
tions with two and three parameters by comparison with 
results of their FEM calculations. As it was found in Ref. 
[J the use of three parameter functions leads to a lower 
error. It was estimated as about 1% when the distance 
from semiconductor/insulator interface is equal to a* B . 

However, as it was shown in Ref. the errors in 
energy values calculated even with the three parameter 
function become much greater in the case Q* = — 1 (semi- 
conductor/metal interface). These errors grow when the 
donor location tends to the interface. The ground en- 
ergy values were calculated in Ref. [l| by FEM, with 
three-parameter and two-parameter trial functions. They 
were determined for R = a* B as equal to Eo/Ry* — 
0.304, 0.284, 0.276, respectively. At the same time the ap- 
proach used in Ref. d gives the value of E /Ry* = 0.3048 
(with the use of 12 functions in our Atrial) which is rather 
close to our FEM value E /Ry* = 0.3050. 

For the system under consideration a simple single- 
parametric function has been suggested earlier in Ref. [4[ ■ 
This function allows to calculate all energies analytically, 
but provide less accurate energy values as compared with 
the results in Ref. [l[ . Another type of a two-parametric 
trial function has been recently suggested in Ref. 
This function allows to perform main calculations analyt- 
ically and gives relative errors for energy values less than 
1.4% for any donor distances from the interface. However 
all these trial functions do not provide the 'spectroscopic 
accuracy' of calculated energies. 

In conclusion, we confirm that the theoretical approach 
suggested by MacMillen and Landman 2 gives rather ac- 
curate energy values for the system of hydrogenic donor 
near the interface between semiconductor and insulator. 
One can use their results as a reference in searching for 
other types of trial functions considering electronic prop- 
erties of such systems. 
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